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ABSTRACT 


An  advanced  numerical  scheme  for  computing  the  free  surface  flow  about  a 
ship  is  presented.  The  present  approach  is  a  higher  order  panel  method  which 
employs  a  parabolic  quadrilateral  as  a  basic  element.  It  can  represent  an  actual 
curved  ship  hull  surface  more  accurately  with  a  smaller  number  of  panels  than  is 
possible  with  a  flat  panel  approximation.  Singularity  strength  can  vary  linearly  or 
quadratically  across  the  panel.  A  sample  unbounded  flow  problem  with  an  analytic 
solution  was  considered  first  to  validate  the  present  approach  and  examine  the 
characteristics  of  accuracy  and  convergence  compared  with  those  of  the 
conventional  low-order  approaches  (constant  source  strength  across  fiat  panel).  The 
Wigley  hull  and  the  Series  60.  Cg=0.60  hull,  which  have  been  studied 
extensively,  were  selected  for  sample  free  surface  flow  computations.  Comparisons 
between  calculated  and  experimental  results  show  good  agreement.  A  modified 
body  boundary  condition  was  suggested  to  simulate  a  dry  transom  stern 
numerically,  and  calculated  results  show  good  correspondence  with  test  data. 


ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  the  General  Hydrodynamics  Research  Program,  administered  by  the  David 
Taylor  Research  Center  (DTRC)  Ship  Hydromechanics  Department,  and  was  funded  by  the  Naval  Sea  Systems 
Command  (NAVSEA),  Hull  Research  and  Technology  Office.  Funding  was  provided  under  Program  Element 
61153N,  Task  Area  SR  0230101,  and  Work  Unit  1522-070.  Additional  funding  was  provided  to  the  third 
author  by  the  National  Science  Foundation  under  Grant  Number  ECS  8515174. 


INTRODUCTION 

Many  attempts  to  predict  the  wave-making  resistance  of  surface  ships  have  been  made  since  Michell's' 

pioneering  work  at  the  end  of  the  19th  century.  However,  due  to  a  variety  of  difficulties,  none  of  the  existing 

methods  have  been  fully  satisfactory.  In  the  present  work  a  comprehensive  modern  effort,  using  a  higher-order 

curved  panel  method,  is  adopted  and  applied  to  the  computation  of  the  free  surface  flow  around  a  surface  ship. 

2 

This  choice  was  partially  motivated  by  a  limitation  in  the  Hess  and  Smith  approach  observed  in  existing 
computational  implementations.  The  Hess  and  Smith  approach  uses  fiat  panels  v;ith  constant  source  strength 
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over  each  panel.  Methods  employing  this  feature  are  quite  sensitive  to  the  geometric  arrangement  of  panels,^ 
Numerical  problems  arise  when  the  solution  is  dependent  on  the  variation  of  panel  size  and  arrangement  for  a 
given  configuration.  Such  problems  place  severe  limitations  on  the  development  of  automatic  panelling  or  other 
complementary  aids,  which  are  important  for  ship  designers  in  that  they  relieve  users  from  the  large  amount  of 
hard  work  and  numerous  judgements  associated  with  producing  the  input  geometry  needed  for  accurate 
numerical  solutions.  Though  the  versatility  of  the  Hess  and  Smith  approach  has  been  appreciated  over  the  past 
decade,  the  need  for  basic  improvement  has  become  increasingly  evident.  In  the  present  work  we  have 
employed  a  curved  panel  as  a  basic  element.  It  is  anticipated  that  the  arrangement  using  curved  panels  will 
require  a  smaller  number  of  curved  panels  than  would  flat  panels  to  approximate  the  actual  surface  of  the  ship 
with  the  same  degree  of  accuracy.  Also,  the  source  strength  distribution  on  the  body  has  been  assumed  to  vary 
linearly  across  each  panel.  For  the  solution  of  the  free  surface  flow  problem  we  distributed  either  sources  or 
doublets  on  the  free  surface.  In  the  doublet  network  we  have  assumed  quadratically  varying  doublet  strengths 
across  each  panel  for  the  sake  of  consistency  with  its  quadratic  geometric  description.  For  the  purpose  of  a 
later  extension  to  the  full  nonlinear  free  surface  boundary  conditions,  the  software  for  curved  free  surface 
panels  has  also  been  included  as  a  runtime  option. 

The  capability  of  the  present  approach  was  carefully  examined  by  applying  it  to  the  computation  of 

unbounded  potential  flow  problems  first  and  then  extending  it  to  free  surface  flo\\  problems.  Comparison  of  the 

numerical  results  of  the  unbounded  potential  flow  around  a  sphere  with  the  exact  solution  showed  that  the 

2 

present  higher-order  curved  panel  method  gives  a  convergence  rate  of  0(h  ).  while  the  rate  of  the  zero-order 
flat  panel  method  is  only  0(h),  where  h  denotes  a  characteristic  length  scale  of  the  panel.  Estimation  of 
numerical  accuracy  of  (he  calculated  results  for  a  free  surface  flow  problem  is  far  more  difficult.  First,  there  is 
no  known  analytical  result  with  which  to  compare.  Other  difficulties  are  more  specific  and  concern  a  variety  of 
numerical  approximations  involved  in  the  treatment  of  boundary  conditions  and  the  infinite  computational  fluid 
domain.  Finally,  one  important  quantity  of  interest,  the  wave-making  resistance,  is  very  small  compared  with 
other  quantities,  such  as  the  velocity  potential  and  the  free  surface  elevation.  In  this  report  emphasis  is  placed 
on  the  accurate  prediction  of  the  free  surface  wave  profile  along  a  ship  hull  instead  of  on  the  wave-making 
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resistance.  The  prediction  of  the  wave  profile  implies  the  detailed  description  of  local  flow  phenomena  while  the 
wave-making  resistance  is  a  global  quantity  usually  obtained  by  integrating  the  pressure,  which  is  a  local 
quantity,  along  the  hull  surface.  It  may  be  fruitless  to  attempt  a  computational  method  for  an  accurate 
prediction  of  the  wave-making  resistance  without  first  achieving  a  corresponding  accuracy  of  the  local 
phenomena.  It  is  also  noted  that  none  of  the  existing  numerical  schemes  have  been  systematically  verified  for 
convergence  and  for  the  numerical  approximations  employed.  The  results  have  been  shown  to  vary  dramatically 
depending  on  such  factors  as  the  panel  density  and  the  numerical  truncation  of  the  computational  fluid 

3 

domain.  In  this  report  the  numerical  results  are  carefully  examined  and  the  results  are  found  to  be  quite 
robust  with  respect  to  changes  in  the  truncated  computational  domain  and  both  the  ship  hull  and  free  surface 
panel  density. 

The  numerical  modeling  of  the  free  surface  flow  around  a  surface  ship  with  a  dry  transom  stern  is  an 
another  major  challenge  for  a  modern  free  surface  problem  code.  One  of  the  outstanding  features  of  the 
present  analysis  scheme  SWIFT(Ship  Wave  Inviscid  Flow  Theory)  is  the  treatment  of  the  transom  stern  hull 
characteristic.  This  feature  is  particularly  important  for  the  Navy  because  most  fast  combatant  ships  have 
transom  stern  hulls.  In  addition,  the  numerous  efforts  for  application  of  SWIFT  to  ship  design  and  per- 
-formance  evaluation  problems  have  motivated  the  development  of  extensive  pre-processor  codes  (PANHULL) 
and  post-processor  software  (PLTPRG)  for  geometry  definition  and  graphic  presentation. 

Validation  computations  for  the  free  surface  flow  problem  have  been  made  for  the  Wigley  hull  and  the 
Series  60  hull,  Cg  =  0.60,  for  which  extensive  experimental  measurements  are  available.  Transom  stern  ship 
hulls  have  also  been  considered  in  detail.  Here,  emphasis  has  been  given  to  the  numerical  results  of  wave 
profiles  and  comparisons  with  the  experimental  measurements.  Predicted  longitudinal  wave  cuts  away  from  the 
ship  have  also  been  compared  with  the  measurements. 
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MATHEMATICAL  FORMULATION 


We  consider  a  ship  which  moves  in  the  positive  x-direction  with  constant  forward  speed  U.  It  is  convenient 
to  introduce  a  Cartesian  coordinate  system'?  =  (x,y.z)  translating  with  the  ship  as  shown  in  Fig.  1.  The  z=0 
plane  is  taken  as  the  undisturbed  free  surface,  the  positive  x-axis  in  the  direction  of  the  ship's  forward  velocity, 
and  the  positive  z-axis  upward.  The  fluid  is  assumed  to  be  ideal,  incompressible  with  constant  density  p .  and 
its  motion  to  be  irrotational.  With  these  assumptions  the  fluid  velocity  vector  V(  T  )  is  represented  as  the 
gradient  of  the  velocity  potential  ♦(  ? )  which  satisfies  the  Laplace  equation . 

7  ^♦(‘x  )  =  0  (1) 

throughout  the  fluid  domain.  On  the  submerged  portion  of  the  ship’s  wetted  surface  S,  the  normal  velocity  of 
the  ship  is  equal  to  that  of  the  adjacent  fluid  such  that 

( T)  =  0  on  S.  (2) 

On  the  free  surface  z  =  C(x,y),  where  C  is  the  free  surface  elevation,  the  kinematic  boundary  condition  is 
expressed  in  terms  of  the  substantial  derivative  D/Dt  =  9/3t  +  V*7  .  in  the  form 

(D/Dt)(  C  -  z  )  =  0  on  z  =  C  •  O) 

An  additional  dynamic  boundary  condition  is  obtained  by  requiring  that  the  pressure  on  the  free  surface  be 
atmospheric; 

(l/2)7^*y^  +  gz  =  0  on  z  =  £,-  (4) 

Since  Eq.  (4)  holds  on  the  free  surface,  its  substantial  derivative  can  be  set  equal  to  zero.  This  derivative  leads 
to  an  exact  nonlinear  free  surface  boundary  condition. 

(1/2)7+  +  g+^  =0  on  z  =  C.  (5) 
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with 


C  =  (l/2g)74i*y^  on  z  = 


(6) 


Finally,  energy  consideration  requires  that  the  velocity  potential  approach  the  uniform  onset  flow  potential  and 
that  there  be  no  waves  far  upstream  of  the  ship: 


where  r  = 


♦  = 


(  -Ux  +  0(1 /r) 


-Ux  +  0(1 /r) 


{x^  + 


2  ^  2.1/2 

y  +  z  } 


X  >  0 
X  <  0 


(7) 


The  set  of  Eqs.  (l)-(7)  formulates  the  exact  boundary-value  problem  within  the  limitations  of  an  ideal 
incompressible  fluid.  However,  the  nonlinear  firee-surface  condition  precludes  any  solution  within  the  present 

4 

linear  framework.  In  the  present  linear  approach  we  followed  Dawson's  approximation  that  the  disturbance 
due  to  free  surface  waves  is  small  compared  to  the  corresponding  double-model  solution.  The  total  velocity 
potential  is  now  expressed  as  the  sum  of  two  potential  functions  t  and  . 


♦  =  ♦  -h  ♦’ 


(8) 


where  4  and  4’  are  the  double-model  velocity  potential  and  the  disturbance  potential  due  to  free  surface  waves, 
respectively.  Under  the  assumption  that  4’  is  small  compared  to  4.  the  free  surface  condition  in  Eq.(5)  can  be 
linearized  with  respect  to  4’  such  that 


[  4  V  ^4-4  ^)  +  4  ^(4  ^-44  +  24  (4  4’  +4  4')  +  24  (4  4'  +4  4’  )  -4 

‘  X  X  y  X  y  X  y  y  x  x  x  y  y  y  x  x  y  y 


4’  (4  ^  +  4  +  4’  (4  ^  +  4  ^)  ]  +  g4'  =0 

X  X  'j  '  y'  X  y  '  ^  z 


C  =  (l/2g)[U^  -4^^  -4y^  -24^4’^  -2*y^'y] 


on  z  =  0 


on  z  =  0. 


(9) 

(10) 


The  simplified  form  of  the  linearized  free  surface  condition  expressed  in  terms  of  double-model  streamline 
coordinates  is  written  as 
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on  z  =  0 


(11) 


+  =2^*11 

C  =  (l/2g)[  ]  on  z  =  0  (12) 

where  I  denotes  the  double-model  streamline  coordinate  and  subscript  I  the  derivative  along  the  l-direction.  The 
mathematical  formulation  then  consists  of  the  Laplace  equation.  Eq.  (I).  with  the  boundary  conditions  of  Eqs. 
(2)  and  (1 1)  on  a  truncated  region  of  the  free  surface  near  the  ship,  and  Eq.  (7).  Equation  (12)  is  used  for  the 
wave  height  computation. 


NUMERICAL  METHOD 

The  method  used  here  is  a  boundary  element  method  in  which  simple  Rankine  sources  are  distributed 
across  each  panel.  The  boundary-value  problem  formulated  above  then  reduces  to  a  determination  of  an 
unknown  singularity  distribution  over  the  boundary  surface  of  the  fluid  domain.  Once  singularity  distribution  is 
determined  as  the  solution,  the  hydrodynamic  quantities  of  interest,  the  velocity  and  the  pressure,  are  also 
determined.  The  integration  of  the  hull  surface  pressure  gives  forces  and  moments  on  the  ship. 

The  new  and  original  features  of  the  present  approach  are  summarized  as  follows: 

1 .  The  geometry  of  all  the  boundary  surfaces  is  patched  by  quadratically  varying  curved  panels  instead 

of  flat  panels. 

2.  Linearly  varying  sources  or  quadratically  varying  doublet  singularities  are  distributed  on  the  surface 

panel,  while  existing  approaches  use  constant  sources  on  each  panel. 

3.  A  new  transom  stern  condition  is  used. 

4.  An  automatic  method  of  determining  a  truncated  free  surface  region  and  a  suitable  density  of  free 

surface  panels  depending  the  Froude  number  is  implemented. 

5.  A  variety  of  graphical  post-processors  have  been  developed. 

The  following  sections  give  a  detailed  description  of  the  advantages  of  the  present  approach  over 
conventional  approaches. 
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HIGHER-ORDER  CURVED  PANEL 


Wc  use  curved  panels  instead  of  flat  ones  to  patch  the  ship  surface.  Let  S  be  a  true  surface  element 
bounded  by  four  corner  points  PI,  P2.  P3.  and  P4  as  shown  in  Fig  2.  We  assume  that  S  may  be  represented 
in  the  approximate  form 

C(t.n)  =  Cjj  +  +  c^^n^/2  (i3) 

where  (^,11,0  arc  orthogonal  coordinates  local  to  S.  The  six  coefficients 

obtained  by  requiring  that  the  approximate  surface  given  by  Eq.  (13)  pass  through  its  four  corner  points  exactly 
and  through  its  neighboring  points  approximately  in  a  least  squares  sense.  That  is.  for  each  panel  the 
expression 


R  =  1/2  1  W.  (H) 

j*l  J  J  J  J 

is  minimized  with  respect  to  the  six  unknown  coefficients  appearing  in  Eq.  (13).  The  summation  in  Eq.  (14) 
ranges  over  the  N  neighboring  grid  points.  The  choice  of  the  N  points  and  the  weights  Wj  has  been  made  in 
an  attempt  to  minimize  any  irregularities  that  may  appear  in  the  paraboloid  approximation  of  the  true  surface. 


In  the  present  formulation  the  values  N  =  9,12,  or  16  and  =  1  or  10  arc  used  depending  on  the  location 
of  the  panel  in  the  surface  panel  network.  A  typical  3-D  perspective  view  of  the  surface  panel  network  is  shown 


in  Fig.  3.  Here  we  use  the  term  network  to  represent  the  collection  of  quadrilaterals.  The  use  of  curved  panels 

reduces  the  geometric  discontinuity  (leakage)  experienced  by  the  usual  faceted  flat  panel  approximations. 

2 

Furthermore,  the  choice  of  these  features  is  required  by  (he  argument  given  by  Hess,  that  certain  computing 


methods  which  use  flat  panels  do  not  obtain  increased  accuracy  over  the  basic  zero-order  method  (constant 
source  on  flat  panels)  even  though  they  use  higher-order  source  distributions.  It  is  shown  by  Johnson^  that  for 
purposes  of  computational  efficiency  Eq.  (13)  can  be  reduced  by  an  appropriate  coordinate  transformation  to 


the  canonical  form 
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C  =  +  bn^ 


(^.M)  E  Z 


(15) 


where  Z  is  the  quadrilateral  formed  by  the  projection  of  the  corner  points  on  the  (£,.f))  plane.  In  the  case  of  a 
flat  panel,  a  and  b  are  both  zero. 

Let 

8  =  1/2  Max  {  a^e^  +  bW  .  (^.n)  e  E  (16) 

and  assume  that 

8  <<  1.0.  (17) 

Physically,  this  assumption  implies  that  the  surface  paneling  is  fine  enough  to  assume  a  small  ratio  of  the 
height  to  the  maximum  diameter  of  the  panel.  Under  this  assumption  the  integrals  associated  with  the  induced 
potentials  and  velocities  are  expressible  as  integrals  in  the  E-plane,  and  all  the  resultant  influence  coefTicients 
are  obtainable  in  closed  form.  Furthermore,  ail  the  influence  coefficients  associated  with  the  higher-order 
curved  panel  method  can  be  computed  from  the  corresponding  integrals  pertinent  to  the  zero-order  flat  panel 
methods.  Since  these  additional  computations  are  carried  out  through  a  linear  recursion  relationship,  the  extra 
effort  is  minimal.  The  justification  of  this  approximation  was  given  by  Johnson.^  His  numerical  tests  suggested 
that  8  <0.()66ensures  a  sufficiently  fine  paneling,  and  that  a  coarser  paneling  with  a  higher  ratio  gives  results 
with  less  accuracy. 

The  distribution  of  singularity  strengths  on  a  panel  S  is  assumed  to  be  linear  in  the  case  of  a  source 
network  and  quadratic  in  the  case  of  a  doublet  network.  Specifically,  we  assume  that  the  singularity  strengths  o 
and  T  at  a  point  on  a  panel  S  are  given  by 

Source:  (18) 

Doublet:  T(£;.n)  “  (>9) 
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The  coefficients  in  Eqs.  (18)  and  (19)  are  determined  by  using  a  method  of  weighted  least  squares  over  the 
panel  and  up  to  eight  of  its  neighboring  panels.  This  method  requires  that  the  form  of  Eq.  (18)  or  (19)  gives 
exactly  the  singularity  value  at  its  centroid  and  approximate  values  in  the  least  square  sense  at  centroids  of  the 
neighboring  panels,  as  shown  in  Fig.  4.  Since  the  singularity  strengths  are  not  known  a  priori,  these 
coefficients  are  obtained  as  linear  combinations  of  the  unknown  singularity  values  at  discrete  centroidal  points. 
For  each  panel  L  there  exists  a  relationship. 


Here  the  subscript  L  identifies  the  panel  of  interest  and  the  subscripts  1,2.  ? . 9  the  local  panel  indices 

associated  with  the  panel  L. 

The  sole  reason  for  this  least  squares  procedure  is  to  express  the  assumed  panel  singularity  distribution  in 
terms  of  a  neighboring  subset  of  the  unknown  singularity,  such  that  this  approach  requires  no  additional 
control  points  to  those  used  in  a  constant-source,  flat-panel  method.  It  still  gives  the  full  power  of  the 
higher-order  singularity  method. 

INDUCED  POTENTIAL  AND  VELOCITY 

Now  the  velocity  potential  at  P“(x,y,z)  induced  by  a  singularity  distribution  on  S  is  given  by 
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Source:  ♦  =  (-l/4n)  JJ  (a/R)  ds  . 


(21) 


Doublet;  4*  =  (-l/4n)  JJ  T(R-n/R^)  ds  . 


(22) 


where 


R  =  (^x,n-y.C-z)  *  0  -T*". 

(23) 

R  =  |R1  =  {  (^x)2-(-0>y)^-l-(C-z)2  . 

(24) 

and  a  and  T  are  given  in  Eqs.  (18)  and  (19),  respectively. 

The  velocity  vector  V  induced  at  P  is  given  as  the  gradient  of  the  corresponding  velocity  potential. 

v"  = 

(25) 

We  are  now  in  a  position  to  evaluate  the  integral  Eqs.  (21).  (22).  and  (25).  From  the  relations  of  Eqs. 
(15),  (17),  (18),  and  (19)  the  induced  velocity  potential  at  P(x.y,z)  due  to  the  singularity  distribution  on  the 
panel  L  can  be  expanded  in  closed  form  by  using  a  combination  of  closed  form  calculations  and  recursive 
relationships.  Details  are  fully  included  in  Johnson’s  report^  and  in  Appendix  A. 

NUMERICAL  APPROXIMATION  OF  THE  FAR-FIELD  BEHAVIOR  OF  THE  SOLUTION 

In  the  present  approach  a  finite  computational  domain  is  used  for  the  numerical  computation  by  truncating 
the  infinite  fluid  domain  at  a  finite  distance  from  the  ship.  If  the  truncation  boundary  is  located  sufficiently  far 
from  the  ship,  so  that  the  contribution  resulting  from  this  truncated  region  is  negligible,  this  numerical 
approximation  can  be  justified.  However,  unless  we  specify  legitimate  boundary  conditions  on  this  far 
boundary,  the  solution  of  the  boundary-value  problem  is  not  unique  and  the  numerical  solution  cannot  be 
justified.  In  order  to  render  the  boundary-value  problem  unique,  or  alternatively  to  ensure  that  the  waves 
associated  with  the  disturbance  appear  only  downstream,  three-point  or  four-point  upstream  difference  schemes 
have  been  commonly  used  for  the  approximation  of  the  terms  t||  and  (^^|^|)|  ^9-  (I  D- 
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(26) 


(3w/3l).  =  C4.C0.  +  C4.  ,0).  ,  +  C4.  -OX  -  +  C4.  -OO.  , 
1  II  1-1  1-1  1-2  1-2  1-3  1-3 

or 

(3h)/31).  =  C3.(aX  +  C3.  ,w.  ,  +  C3.  _ 

II  1-1  1-1  1-2  1-2 


(27) 


where  the  coefficients  C4  and  C3  are  given  by 


where 


or 


where 


and 


C"i-1 

=  -‘‘'i-2‘*'i-3<‘“i-3-‘‘'i-2>/°" 

(28) 

C4i-2 

=  ‘»i-l‘»i-3(‘»i-3-‘‘'i-P'®" 

(29) 

C^i-3 

=  -‘»i-l‘"i-2^‘“i-2-‘*'i-l>'^" 

(30) 

C4. 

1 

=  -(C4..,  4-  C4..2  +  C4..3) 

(31) 

D4  = 

‘•'l-l‘''i-2‘‘h-3<‘‘V3-"Vl><‘>h-3*‘*'i-2><‘*‘i-2*‘''^^ 

(32) 

C3i  =  (dl.,2-dl..j)/D3 

(33) 

1 

u 

=  -dl.  -/D3 

1-2 

(34) 

^^1-2 

=  dl.  ,/D3 

1-1 

(35) 

D3  «  dl.  ,dl.  -(dl.  ,+dl.  -) 

i'1  1-2  1-1  1-2 


dl.  ,  -  I.  -  I.  , . 

i-k  I  i-k 


(36) 

(37) 


This  upstream  finite  difference  scheme  for  the  frec-surface  boundary  condition  introduces  numerical 
dissipation  into  the  physical  problem  to  ensure  only  downstream  waves.  However,  the  extent  of  this  numerical 


dissipation  depends  on  the  panel  sizes  and  on  the  free-surface  panel  arrangement  so  that  it  is  hard  to  control  or 
to  estimate.  Furthermore  the  upstream  difference  scheme  introduces  unnecessary  numerical  dispersion  into  the 
physical  problem.  An  alternative  way  to  ensure  no  upstream  waves  is  to  introduce  explicitly  an  artiricial 
viscosity,  known  as  Rayleigh  viscosity,  so  that  the  numerical  dispersion  due  to  the  difference  scheme  can  be 
avoided.  The  corresponding  free-surface  boundary  condition  becomes 


[  ♦  +  2#  ♦  tj;  +  ♦  2(i  ♦  +♦  ♦  )\)/  +  2{t  ♦  ♦  )t|/  ] 

X  ^xx  X  y’^xy  y  yy  x  xx  y  xy  x  x  xy  y  yy  y 


-I-  V  [(2*^-1-  ]  +  gtk^ 

•=  2t  ♦  -t-  2*  ♦  ♦  +  -h  *  *  -ft*  ]  +  V  [(*  +1)^  -I-  #  1  (38) 

‘xxx  xyxy  yyy  yxy  x  y  ^  ' 


where  V  is  a  small  positive  Rayleigh  viscosity  term  and  tk  is  the  disturbance  potential  such  that  *  -Ux  -t-  t|;. 
The  derivative  terms  in  Eq.  (38)  are  evaluated  analytically  through  the  recursion  relationship  with  minor 
additional  computations.  The  dissipation  introduced  by  an  artificial  Rayleigh  viscosity  should  be  small  enough 
not  to  damp  the  solution  too  much. 

The  radiation  condition  of  no  waves  upstream  is  imposed  to  make  the  solution  unique.  If  we  had  applied 
appropriate  boundary  conditions  on  this  far-fleld  truncation  boundary  to  ensure  uniqueness,  the  approximate 
procedures  could  have  been  avoided.  A  physically  acceptable  uniqueness  condition  is  the  condition  of  no  waves 
at  the  upstream  truncation  boundary.  That  is. 


'"x 


(39) 


♦  -  0. 

XX 


(40) 


Once  the  singularity  distributions  are  determined  as  the  solution  of  the  boundary-value  problem  formulated 
above,  the  hydrodynamic  pressure  coefficient  can  be  obtained  from  the  Bernoulli  equation 


Cp  -  p/(l/2  pU^)  -  1  -  V^/U^. 


(41) 


A  linear  distribution  of  C  on  each  panel  can  be  calculated  from  the  values  of  C  at  the  control  points.  The 

P  P 
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method  used  is  identical  to  that  for  computing  the  linear  source  distribution  network.  On  any  panel  S.  we  then 
have 

c  (t.n)  =  c  +  c  +  c  (42) 

p"^’  "  po  pt  pl1 
The  force  coefficient  vector  Cp  on  S  is  defined  by 

Tp  =  (1/S)  JJ  Cp(t,n)n(tlT)ds  (43) 

where  S  is  a  specified  reference  area.  After  algebraic  manipulations,  we  obtain  the  following  expression  for  Cp 
in  local  coordinates: 

^P  =  (1/S)  {  -2al  C(2,l)Cpp  +  C(3.1)Cp^  +  C(2,2)Cp^  ]. 

-2b[  C(1.2)Cp^  +  C(2,2)Cp^  +  C(1.3)Cp^]. 

[  Cd.DCp^  +  C(2.1)Cp^  + 

IMPLEMENTATION  OF  NUMERICAL  SOLUTION 

The  present  method  uses  a  discrete  panel  method  to  achieve  a  numerical  solution  of  the  boundary-value 
problem.  More  specifically,  the  problem  is  reduced  to  the  determination  of  unknown  singularity  strengths  at 
discrete  centroidal  points  on  each  panel.  A  finite  set  of  control  points  (equal  in  number  to  the  number  of 
unknown  singularity  strengths)  is  selected  at  which  the  boundary  conditions  arc  imposed.  Application  of  the 
boundary  conditions  results  in  a  finite  set  of  simultaneous  algebraic  equations  denoted  symbolically  by 

[A]  {X}  -  {B}.  (45) 

Here,  {X}  Is  the  vector  of  unknown  singularities.  {B}  is  the  vector  of  specified  boundary  conditions,  and  [A] 
is  the  square  matrix  of  influence  coefficients.  A.j  represents  the  influence  of  the  j-th  singularity  X.  on  the 
boundary  condition  at  the  i-th  control  point. 
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While  the  solution  of  Eq.  (45)  may  be  obtained  by  using  any  available  numerical  algorithm,  we  have  paid 
particular  attention  to  the  techniques  of  numerical  linear  algebra  in  this  research.  The  part  of  the  matrix 
resulting  from  the  ship  hull  surface  boundary  condition  is  diagonally  dominant,  which  has  led  many 
researchers  to  use  iterative  methods  analogous  to  aerodynamic  computations.  When  the  linearized  free  surface 
conditions  are  added  for  the  study  of  free  surface  flows,  this  diagonal  dominance  is  lost.  It  is  common  for 
research  purposes  to  solve  free  surface  flows  at  a  number  of  different  ship  velocities  without  altering  the 
paneling  with  fixed  sinkage  and  trim.  Letting  the  number  of  body  surface  panels,  free  surface  panels,  and  the 
total  panels  be  denoted  by  NBP,  NFP,  and  NTOT,  it  is  easy  to  sec  that  the  resulting  matrix  is  of  size 

(NTOT)^  =  (NBP+NFP)^.  (46) 

This  matrix  may  be  partitioned  as 


'^B 

'Sf 

ApB 

ApF 

(47) 


Only  the  submatrices  Ap^  and  App  change  with  ship  speed.  A^^  and  App  being  generated  from  the 
no-penetration  condition  on  the  wetted  hull  surface.  A  partitioned  matrix  solution  technique  can  be  deviced  to 
eliminate  a  large  number  of  redundant  calculations,  requiring  matrix  work  onh  on  App  and  App.  A  major 
benefit  of  this  partitioned  matrix  approach  is  that  only  the  submatrices  need  to  be  in  core  at  any  given  stage  of 
computation,  which  greatly  increases  the  number  of  ship  surface  panels  and  especially  the  number  of  free 
surface  panels  available  to  the  user.  All  the  experience  of  the  aerodynamic  researchers  has  indicated  an 
increasing  demand  for  more  panels  to  enable  the  study  of  more  complicated  geometries.  In  naval  applications 
also  a  large  number  of  free  surface  panels  is  often  desired  and  a  total  of  mound  1000  can  be  a  serious 
limitation.  With  the  program  developed  here,  the  only  limitations  on  size  are  that  number  of  the  body  surface 
panels  and  the  free  surface  panels  must  each  individually  be  not  more  than  900  for  use  on  a  CRAY  X'MP/12. 
and  this  limit  increases  to  over  1400  on  a  CRAY  X-MP/24  such  as  >s  available  at  DTRC.  Thus  the  total 
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number  of  panels  will  be  between  1800  and  2800  on  the  CRAY  X-MP  family  of  supercomputers,  depending  on 
available  memory  sources. 

Partitioning  the  solution  vector  {X}  and  the  right-hand  side  vector  {B}  in  a  compatible  way  gives  the  total 
partitioned  matrix  equation 


^B  "^BF 

("M  -  1 

(Bb| 

^FB  "^FF  i 

Bp 

. 

This  equation  reduces  to 

f  App  -  Apg  fAgg  Agp]  ]  {Xp}  =  {Bp}  -  Apg  [  Agg  Bg] 


(48) 


(49) 


with  Xg  being  g'ven  by  an  almost  no-cost  calculation  as 

{Xg}  =  (Agg  Bg)  -  [  Agg  Agp  ]  {Xp}. 


(50) 


Note  that  Agg**  Agp  and  Agg'*Bg  can  be  computed  just  once  so  that  the  only  new  work  of  consequence 

needed  for  additional  ship  speed  cases  is  the  above  matrix  times  matrix  multiplication  and  a  Motorization  of  the 

modiPied  NFP  by  the  NFP  free  surface  matrix  in  Eq.  (49).  If  the  solution  for  only  one  Froude  number  is 

3 

required,  the  total  operation  count  is  the  usual  (NTOT)  /3  exactly  as  when  no  partitioning  is  used,  except  that 
the  operations  for  the  Agg  factorization  are  not  repeated. 

We  have  also  implemented  and  extended  some  very  recent  advances  in  squeezing  supervector  speed  out  of 
matrix  multiplication  and  factorization  routines  recently  developed  by  Dongarra  and  Eisenstat.^  This  procedure 
improved  the  vectorization  in  the  present  numerical  code  so  much  that  the  resulting  factorization  times  in 
portable  FORTRAN  were  signiFicantly  better  than  the  previously  used  LINPACK-based  code  optimized  in 
assembly  language. 

Numerous  other  vectorization  techniques  were  used  which  were  variations  on  the  theme  of  creating  large 
inner  loops.  In  some  cases  reorganizing  algorithms  actually  decreased  the  number  of  operations  significantly. 
For  example,  a  combination  of  both  techniques  decreased  the  computer  time  spent  in  the  computation  of  the 
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far-field  approximation  of  influence  coefficients  by  a  factor  of  85.  An  example  of  the  first  technique  arises  in 
the  typical  computation 

TRJOj^  =  TRJOj^  +  C*TRJj,  .  I  <  K  <  N  (31) 

where  C  is  a  fixed  3  by  3  coordinate  transformation  matrix  and  TRJ^  and  TRJOj,  are  3  by  3  matrices.  This 
calculation  is  routinely  performed  in  the  middle  of  a  large  K-loop  by  calling  on  a  utility  routine  MXM  to 
perform  the  3  by  3  matrix  multiplication  of  C  times  TRJj,.  The  code  can  be  reorganized  so  that  the  K-loop  is 
around  just  this  calculation,  and  a  vector  subroutine  can  be  created  with  the  result  given  by  nine  hand-coded 
dot  products  of  vectors  of  length  3.  One  of  the  nine  FORTRAN  lines  would  read 

TRJO(  1 , 1  ,K)  =  TRJO(  1 , 1  ,K)  -f  C(  1 . 1  )*TRJ(1 . 1  .K)  -I-  C(  1 .2)*TRJ(2 , 1  .K )  -I-  C(  1 .3)*TRJ(3 . 1  .K) . 

This  new  routine  with  the  key  calculation  inside  a  large  loop  is  so  efficient  that  it  decreases  the  computer  time 
for  vectors  of  length  K  around  250  or  larger  by  a  factor  of  15.  The  current  code  is  now  well  structured  for 
inclusion  of  new  modeling  approaches  and  is  over  90  percent  vectorized  for  extremely  efficient  use  on  Class  VI 
and  Class  VII  supercomputers.  It  consists  of  about  12.000  lines  of  FORTRAN.  40  percent  of  which  arc 
comments. 

TRANSOM  STERN 

Most  Navy  high  speed  combatant  ships  have  transom  sterns  and  at  high  speed  the  flow  clears  the  stern 
completely.  To  model  this  unwetted  transom  flow,  one  suggested  approach  is  to  add  a  cruiser  stern-like 
extension  behind  the  transom  stern  and  apply  the  boundary  condition  =  0  on  the  transom  section.  This 
approach  has  given  results  that  compare  poorly  with  experimental  data.  Adequate  numerical  modeling  of  the 
transom  stern  condition  is  challenging  within  the  framework  of  the  linearized  free  surface  potential  flow 
theory.  In  the  basic  formulation  assumptions  have  been  made  that  the  disturbance  created  is  small  enough  to 
allow  the  linearization  of  the  free  surface  condition  about  the  undisturbed  free  surface  C  =0.  Several  variations 
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of  numerically  modeling  have  been  attempted.  The  version  currently  in  the  program  provides  the  most 
reasonable  results  compared  with  experimental  vahies.  In  this  approach,  we  removed  the  transom  stern  part 
from  the  hull  geometry  input  data,  creating  is  a  hole  at  the  stern  of  a  ship  in  the  sense  that  the  hull  is  not 
completely  closed.  Along  the  first  row  of  bottom  panels  adjacent  to  the  transom  stern,  the  hull  surface 
boundary  condition  =  0  was  replaced  by  the  requirement  that  the  flow  pass  by  each  panel  tangentially.  The 
magnitude  and  angle  of  tangential  velocity  were  obtained  by  extrapolating  values  from  neighboring  panels.  The 
free  surface  behind  the  transom  exists  at  z  =  0  so  that  the  linearized  free  surface  condition  is  satisfied.  This 
numerical  model  has  been  tested  on  the  Athena  hull  form  and  on  several  combatant  ships. 

CONVERGENCE  TEST  OF  NUMERICAL  SOLUTION 

In  addition  to  the  difficulty  experienced  in  the  numerical  implementation  of  the  radiation  condition,  there 
exist  numerous  sources  of  computational  uncertainties,  such  as  (1)  numerical  differentiation  and  integration, 
(2)  truncated  domain  from  infinite  to  finite,  and  (3)  panel  size.  Furthermore  these  sources  of  errors  are 
interwoven  in  a  compilicated  manner.  It  is  not  possible  to  separate  and  quantify  simply  each  source  of  error. 
Three  main  attributes  generate  traceable  errors:  (I)  fore-  and  aft-  truncation  boundary.  (2)  side  truncation 
boundary,  and  (3)  number  of  panels  for  a  given  finite  domain.  These  attributes  are  examined  separately  below. 
A  Wigley  parabolic  hullform  in  Fig.  5  was  used  for  this  numerical  convergence  experiment. 

Fore  and  Aft  Truncation 

Fore  and  aft  truncation  boundaries  are  always  taken  the  same  distance  from  FPfforward  perpendicular)  and 

APfaft  perpendicular)  of  the  ship,  respectively.  Their  ranges  arc  approximately  between  1 .5L  and  3.2L  where 

L  is  the  half  of  the  ship  length,  L*=LWL/2.  Froude  numbers  considered  were  between  0.2  and  0.375.  The 

side  boundary  was  fixed  at  y  =  I  .OL.  The  relative  error  is  defined  as 

[C^  (at  x  =  x.)  -  (at  x  =  .V2L)] 

Relative  Error  =  - 

C  (at  x=3.2L) 
w 

where  C  is  the  wave  resistance  coefficient  and  x.  =  I  .46L.  1 .9L.  2.2L.  2.9L.  and  3.2L.  The  results  arc 
w  I 
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plotted  in  Fig.  6.  The  maximum  relative  error  is  about  2%  and  decreases  proportionally  with  increasing 
distance  of  the  truncation  boundary.  For  this  Fronde  number  range,  the  truncation  boundary  x=2.5L  seems 
sufficiently  far  away  to  provide  stable  results. 

Side  Truncation 

The  examined  side  truncation  boundaries  were  located  between  0.5L  and  1.6L  from  the  ship  s  centerline. 

Fore-  and  aft-  boundaries  were  fixed  at  x=2.0L  and  the  same  Froude  number  range  was  considered.  The 

relative  error  is  defined  similarly  as 

[C  (aty=y)-C  (aty=1.6L)] 

Relative  Error  =  -  - -  —  -  -  ■ 

(at  y=1.6L) 

where  y.  =  0.53L.  l.OL,  1.25L.  and  I.6L.  Figure  7  shows  the  results.  The  observed  maximum  error  is  about 
2.5%  and  decreases  linearly  with  the  distance  of  the  boundary.  The  differences  are  hardly  distinguishable  when 
y  is  more  than  1 .25L. 

Panel  Size 

The  last  parameter  considered  is  the  panel  size  or  number  of  panels  for  a  fixed  domain.  The  computational 

domain  was  fixed  at  x=  2.0L  fore  and  aft  and  y=1.0L.  The  total  number  of  panels  for  both  body  and  free 

surface  varied  from  612  to  1052.  Figure  8  shows  the  relative  error  which  is  the  ratio  between  the  difference  of 

at  612,  700,  800,  860,  950,  and  1052  panels  and  of  1052  panels.  The  maximum  error  is  about  4%. 

The  error  decreases  rapidly  to  about  1%  between  600  and  700  panels,  and  at  more  than  700  panels  the  rate  of 

decrement  is  very  slow  compared  with  the  rate  of  increment  in  panel  numbers.  It  appears  that,  once  the 

primary  variation  is  obtained,  the  effect  of  further  refinement  of  the  panel  size  is  not  significant.  In  problems  of 

this  type,  with  abrupt  changes  in  the  boundary  conditions  (hull-free  surfoce  interface)  or  geometry,  there  is  a 

singularity  in  the  solution  within  a  local  neighborhood.  This  singularity  leads  to  the  slow  convergence  of  the 

solution.  Typically  a  halving  of  the  panel  size  (which  increases  the  work  by  a  factor  of  4)  might  produce  only  a 

* 

40%  improvement  in  the  error.  In  recent  work.  Lucas  and  Oh  formulated  a  ne^v  approach  to  this  problem 

*  "The  Method  of  Auxiliary  Mapping  and  the  Finite  Element  Solution  of  Elliptic  Boundary  Value  Problem 
Containing  Boundary  or  Corner  Singularities  submitted  to  IMA  J.  of  Numerical  Analysis. 
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which  has  promise  for  a  generalization  to  the  3- dimensional  problem.  For  the  moment,  however,  aliout  all  that 
can  be  said  i.^  that  a  medium-sized  grid  give.s  sufficient  resolution.  Without  a  new  and  different  approach, 
further  refinement  is  probably  too  costly  to  be  effective. 

RESULTS  AND  DISCUSSION 

The  first  step  of  the  panel  method  approach  is  to  discretize  the  fluid  boundary  in  the  form  of  panel 
networks.  These  networks  are  a  collection  of  either  source  or  doublet  panels  and  are  independently  defined 
over  various  portions  of  the  ship  hull  and  the  free  surface  as  illustrated  in  Fig.  9.  Once  this  panel  dis- 
-cretization  has  been  completed,  the  solution  of  the  boundary-value  problem  can  be  carried  out. 

SPHERE 

The  first  test  case  is  an  unbounded  potential  flow  about  a  sphere  in  a  uniform  stream.  The  solution  is  well 
known  analytically.  This  choice  was  made  to  examine  the  capability  of  the  present  approach  versus  that  of  the 
zero-order  flat-panel  approach.  There  is  symmetry  of  the  flow  about  two  planes  so  that  only  a  quarter  of  the 
sphere  needs  to  be  paneled.  To  test  the  flexibility  of  the  panel  method  we  used  two  different  panel  arrangements 
as  shown  in  Fig.  10:  a  regular  paneling  and  a  random  paneling.  The  corner  points  for  random  panels  were 
chosen  using  a  random  number  generator,  which  led  to  a  wide  variation  in  panel  size  and  shape:  the  corner 
points  of  regular  panels  were  generated  using  cosine  spacing  along  the  x-direction  and  uniform  spacing  along 
the  transverse  direction.  Velocity  potentials  at  each  of  the  12x4  control  points  are  plotted  in  Fig.  10  as  a 
function  of  the  x-coordinate.  Agreement  with  the  exact  solution  is  good.  This  example  demonstrates  the  extreme 
forgiveness  of  the  method  to  irregular  paneling,  a  feature  which  greatly  enhances  the  method's  practical  utility 
for  applications  involving  much  more  complex  configurations.  The  accuracy  and  the  convergence  rate  of  the 
method  are  also  compared  with  those  of  the  zero-order  flat  panel  approach.  The  results,  summarized  in  Table 
1 ,  show  the  superiority  of  the  present  higher-order  cur\’ed  panel  method  over  the  conventional  zero-order  flat 
panel  method  in  both  accuracy  and  convergence  rate:  the  present  method  gives  a  conN-ergence  rate  of  O(h^). 
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while  the  rate  of  the  conventional  method  is  only  0(h).  Here  h  denotes  a  characteristic  length  scale  of  the 
panel. 

The  results  in  Table  1  for  the  zero-order  flat  panel  approach  were  obtained  as  the  special  limiting 

solutions  of  the  present  higher-order  curved  panel  code.  Numerical  tests  showed  that  the  zero-order  flat  panel 

2 

solutions  from  the  present  code  are  exactly  same  as  those  from  an  existing  zero-order  code. 

WIGLEY  HULL 

The  Wigley  hull  is  a  mathematical  hull  with  a  parabolic  waterline  and  parabolic  cross  section.  The 
equation  used  for  the  surface  is 

y  =  (l-x^)(l-z^/64)/IO  ,  (52) 

and  the  particulars  and  body  plans  are  given  in  Fig.  5.  Six  Froude  numbers.  0.25,  0.267.  0.289.  0.316, 
0.354,  and  0.408,  were  considered.  Extensive  numerical  tests  revealed  that  the  extent  of  the  paneled  part  of  the 
free  surface  must  be  determined  on  the  basis  of  the  wavelength  of  the  generated  waves  as  well  as  on  the  hull 
configuration .  Numerical  tests  suggest  that  a  finite  free  surface  region,  extending  one  wavelength  upstream 
from  the  bow,  one  wavelength  downstream  from  the  stern,  and  3/4  ship  length  from  the  centerline  to  the  side, 
is  sufficient  to  approximate  an  unbounded  free  surflice  to  the  desired  accuracy. 

It  is  a  common  practice  in  the  field  of  computational  fluid  dynamics  to  verify  the  accuracy  of  a  numerical 
solution  by  systematically  varying  panel  densities  and  then  examining  the  resulting  convergence  characteristics. 
Systematic  variations  of  panel  densities  on  the  wetted  hull  surflice  and  on  the  free  surflice  were  investigated  for 
this  mathematical  hull  form.  During  these  numerical  tests  the  upstream  finite  difference  operator  given  in  Eq. 
(26)  or  (27)  was  used  to  approximate  the  radiation  condition  numerically.  The  effect  of  the  free  surface  panel 
density  on  the  numerical  solution  was  then  tested  by  changing  the  densities  in  the  lateral  and  longitudinal 
directions.  It  was  found  that  the  numerical  results  were  more  sensitive  to  the  variation  of  the  free  surface  panel 
density  than  to  the  body  surface  panel  density.  This  sensitivity  may  be  explained  in  two  ways.  First,  the  free 
surface  boundary  condition  is  of  mixed  type  Involving  tangential  velocities  as  well  as  normal  velocities,  while 
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the  body  boundary  condition  is  of  Neumann  type  involving  only  velocities  normal  to  the  boundary.  Secondly, 
we  use  the  upstream  finite  difference  operator  to  compute  the  streamwise  derixatives  of  the  tangential  velocities 
in  the  free  surface  condition.  This  operator  is  used  to  eliminate  upstream  disturbances  numerically.  It  is  well 
known  that  the  accuracy  of  the  finite  difference  operation  depends  on  the  size  of  the  basic  element.  Furtherm¬ 
ore,  the  artificial  numerical  damping  involved  in  a  specific  upstream  finite  difference  operation  is  expected  to 
vary  with  the  free  surface  panel  densities. 

In  the  present  report,  we  used  for  the  higher-order  panel  approach  a  linearly  varying  source  distribution  on 
the  curved  ship  hull  and  a  quadratic  doublet  distribution  on  the  free  surface.  In  order  to  satisfy  the  free 
surface  condition  a  three-point  upstream  finite-difference  and  a  four-point  upstream  finite  difference  operator 
were  used  for  the  higher-order  panel  approach  and  the  zero-order  panel  approach,  respectively.  The  ship  hull 
surface  was  defined  using  192  panels  (24x8),  and  624  panels  (52x12)  w'ere  used  to  represent  the  truncated  free 
surface  domain.  In  Figs.  11  through  17,  the  computed  wave  profiles  and  wave-making  resistance  results,  are 
compared  with  the  measurements  reported  by  the  Resistance  Committee  of  the  1 7th  ITTC(International  Towing 

7 

Tank  Conference). 

To  focilitate  comparisons,  the  following  non-dimensionalization  has  been  made: 

2 

C^(Wave-making  resistance  coefficient)  =  Resistance/(l/2  pU  S) 

C(Wave  profile)  =  0(U^/2g) 

where  U  is  the  forward  speed,  g  the  gravitational  acceleration,  and  S  the  wetted  surface  area,  (generally 
speaking,  the  advantage  of  the  higher-order  panel  approach  is  not  evident  in  these  numerical  results  to  the 
same  extent  that  the  results  for  the  unbounded  potential  flow  show  the  advantage  of  the  higher-order  panel 
approach  over  the  zero-order  panel  approach.  The  numerical  results  shown  in  the  Figs.  1 1  through  17  suggest 
that  a  fair  comparison  of  approaches  must  be  made  w'ithout  using  an  approximate  scheme  (upstream-finite 
difference  operator)  which  induces  numerical  errors  that  cannot  be  estimated  easily.  The  alternative  approaches 
using  Eqs.  (38)  or  (39)  and  (40)  instead  of  the  upstream-finite  difference  operator  are  expected  to  provide  fair 
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comparisons  and  to  show  the  superiority  of  the  higher-order  approach  even  for  the  free-siirface  flow  problems. 
Figure  18  shows  3-dimensional  perspective  view'  of  waves  generated  by  Wigley  hull  at  =0.408.  Wave 
contour  plot  is  shown  in  Fig.  19. 

SERIES  60.  Cg  =  0.60 

The  body  plan  and  profile  drawings  of  the  cruiser  stern  Series  60.  Cg  =  0.60  hull  form  are  shown  in  Fig. 
20  and  the  table  of  offsets  is  provided  in  Table  2.  In  this  case,  the  ship  hull  surface  was  defined  using  360 
panels  (30x12),  and  696  panels  (58x12)  were  used  for  the  truncated  free  surface.  Figures  21  through  26  show 
the  measured  and  the  computed  wave  profiles  compared  at  six  Froude  numbers:  0.20.  0.25.  0.30.  0.32.  and 
0.35.  At  the  forward  half  of  hull  the  computed  wave  profiles  compare  favorably  with  the  measured  ones,  but 
the  agreement  becomes  poorer  downstream.  The  prediction  always  overestimates  the  magnitude  of  the  last 
crest.  The  discrepancy  near  the  stern  may  be  explained  partially  by  the  fact  that  growth  of  a  thick  boundary 
layer  violates  the  underlying  assumption  of  potential  flow,  and  the  nonlinear  free  surface  condition  should  be 
taken  into  account  in  order  to  make  improvement.  Note  that  the  comparison  between  the  computation  and 
experiment  for  the  Series  60  hull  is  different  from  the  Wigley  hull  case  in  which  the  prediction  always 
underestimates  the  magnitude  of  the  bow  crest  and  shows  good  agreement  in  the  aft  part. 

In  Fig.  27,  the  calculated  wave-making  resistance  coefficients  are  compared  with  the  experimental  results 

g 

of  Kim  and  Jenkins.  It  is  interesting  to  note  that  for  both  the  Wigley  and  Series  60  hulls  the  differences 
between  the  measured  and  the  computed  wave  profiles  are  relatively  small,  but  the  discrepancy  in  wave-making 
resistance  between  theory  and  experiment  is  much  greater  than  expected.  As  is  well  known,  most  of  the 
wave-making  resistance  comes  from  the  difference  between  the  integrated  longitudinal  pressure  forces  at  the 
bow  and  at  the  stern.  Particularly  in  these  two  areas,  the  observed  flow  seems  to  violate  the  underlying 
assumptions  for  linearized  potential-flow  theory,  in  Fig.  28  the  computed  longitudinal  wave  cut  was  compared 
with  the  measured  one  at  y/(L/2)“0.752l  and  F^  =  0.323.  Near  and  aft  of  the  stern  the  wave  phase  is  shifted 
slightly,  but  the  amplitude  shows  excellent  agreement  between  calculated  and  experimental  results. 
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ATHENA  HULL 


A  body  plan  of  the  round  bilge,  semi-planing  high  speed  ship  hull  of  the  Athena  is  shown  in  Figs,  29  and 
30.  and  the  offsets  are  given  in  Table  3.  The  hull  was  represented  by  23  by  8  panels.  The  transom  stern  and 
free  surface  region  behind  the  transom  stern  were  considered  as  a  special  section  in  the  computation.  A  total  of 
489  panels  was  used  to  cover  the  whole  free  surface;  48x8  panels  for  the  main  |X)rtion  of  the  free  surface  and 
15x7  panels  for  the  free  surface  behind  the  transom  stern.  As  mentioned  previously  the  transom  stern  segment 
of  the  body  geometry  is  left  unpaneled  and  a  tangential  condition  was  satisfied  at  the  last  row  of  panels  of  the 
hull.  Four  Froude  numbers,  0.350,  0.412,  0.484.  and  0.651.  were  selected  from  among  the  test  runs  of  the 
ITTC  Resistance  Committee,'^  and  the  measured  wave  profiles  were  compared  with  the  computed  ones  in  Figs. 
31  -  34.  At  lower  Froude  numbers  a  discrepancy  is  noticable.  especially  aft  of  midship.  The  computational 
results  show  a  deeper  trough  than  do  the  experimental  data.  Furthermore  up  to  Froude  number  0.412  we  see  a 
kink  in  the  wave  profile  very  near  the  stern  which  does  not  exist  in  the  experimental  results.  This  peculiar 
result  at  lower  Froude  numbers  may  come  from  the  transom  boundary  condition  incorporated  here,  since  at 
lower  speeds  the  transom  stern  is  not  completely  dry.  As  speed  increases,  this  phenomenon  disappears 
completely  and  the  predicted  wave  profiles  smoothly  follow  the  measured  ones  along  the  entire  hull.  These 
numerical  results  indicate  that  the  suggested  transom  stern  boundary  condition  is  well  suited  to  higher  Froude 
numbers  (i.e.,  dry  transom  stern),  but  for  lower  Froude  numbers  we  will  have  to  develop  a  wetted  transom 
version  of  the  stern  boundary  condition. 

Figure  35  shows  the  wave  profiles  behind  the  stern  of  the  ship  at  two  different  locations;  y/(L/2)  *  0.027 
and  0.054  for  F^>b0.484.  The  computed  and  measured  results  show  excellent  agreement. 

LARGE  COMBATANT  SHIPS 

As  part  of  an  intensive  program  to  study  the  wave  svstems  generated  bv  a  large  combatant  ship.  DTRC 

9 

conducted  model  tests  for  three  similar  hull  configurations:  Models  X.  A.  and  U  (see  Beaver's  report  ).  The 
purpose  of  the  testing  was  to  study  the  effects  of  hull  shape  changes  on  the  wave  generation  characteristics. 
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Figure  36  shows  a  comparison  of  body  plans  for  Model  X  and  Model  A.  Model  X  is  the  existing  configuration 
and  Model  A  identifies  the  notched-hull  configuration  in  which  the  7.5  fooi  wide  notch  extends  the  length 
between  stations  3.7  and  18.8  at  the  designed  waterline.  The  body  plan  of  Model  C  is  compared  with  that  of 
Model  A  in  Fig.  37.  Model  C  identifies  a  modification  of  Model  A  that  incorporates  only  the  removal  of 
volume  from  the  shoulder  region  between  stations  4  and  10.  Hull  form  characteristics  of  Models  X.  A.  and  C 
are  given  in  Table  4. 

The  wave  profiles  along  the  hulls  were  measured  at  F^  =  0.179.  0.238.  and  0.298  between  stations  0  and 
8  for  each  configuration.  In  Figs.  38.  39,  and  40  the  computed  and  measured  wave  profiles  are  plotted 
together  to  make  comparison  easier  between  Models  X  and  A  at  each  speed.  In  Figs.  41.  42.  and  43  the  same 
comparisons  are  shown  for  Models  A  and  C.  As  shown  in  these  figures,  the  computed  results  agree  well  with 
the  measured  ones  for  all  three  models  and  for  all  speeds  considered,  except  near  the  bow  region. 

Off-body  wave  elevations  for  Models  X  and  A  were  measured  at  F^  =  0.238  using  five  capacitance  wire 
wave  probes  mounted  transversely  across  the  basin.  As  shown  in  Fig.  44.  the  probe  locations  were  approxi¬ 
mately  77,  90.  102,  120  .  and  135  ft  ship  scale  off  the  centerline  of  the  ship.  The  comparisons  between  the 
computed  and  measured  off-body  longitudinal  wave  elevations  at  each  probe  were  plotted  in  Figs.  45  through 
49  and  their  agreements  are  excellent. 

COMPARISON  WITH  OTHER  APPROACHES 

10 

DTRC  organized  a  special  workshop  called  "Wake-Off"  which  was  held  on  12  and  13  January  1988. 
The  objective  of  this  workshop  was  to  gather  all  available  numerical  codes  that  may  be  useful  for  predicting  or 
evaluating  ship  generated  waves,  wave  spectra,  wave  contour,  wave  pattern,  and  wave-making  resistance  and  to 
make  comparative  evaluation  of  results  of  these  codes  for  given  specific  ships  and  speeds.  Two  models.  DDG5I 
and  OUAPAW.  were  selected  and  their  offsets,  particulars,  and  siiecific  Froude  numbers  were  provided  to  each 
participant.  DTRC  conducted  model  tests  at  the  same  time,  but  the  test  dam  were  not  disclosed  to  the 
participants.  This  closure  provided  a  blind  test  for  predictive  capability  of  the  various  computational 
approaches.  Twelve  computer  programs  were  involved  in  this  workshop:  five  from  DTRC.  two  from 
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SAIC(Scicnce  Application  International  Corporation),  one  each  from  Uni'ersit}'  of  Michigan.  Lockheed.  SSPA 
Maritime  Consulting.  MTT(Massachii setts  Institute  of  Technology),  and  Rosen(Lngineering  Consultant).  It  was 
interesting  to  see  that  only  two  methods,  the  Havelock  source  and  the  Rankine  source  method,  were  used 
among  them.  Five  computer  codes  use  the  Havelock  source  and  the  others  use  the  Rankine  source  method. 

Both  methods  consistently  underpredicted  the  amplitude  of  the  bow  wave  and  o\'erpredicted  the  amplitude  of 
the  stern  wave.  Overall,  the  Havelock  source  codes  provided  better  results  than  the  Rankine  source  codes, 
especially  for  the  far  field  wake.  Excellent  near-fleld  calculations  were  provided  by  the  codes  SWIFT  and 
FLOPAN  which  use  a  higher  order  Rankine  source  and  a  constant  source/doublet  distribution  method, 
respectively.  These  two  codes  outperformed  the  methods  such  as  XYZFS  which  used  the  constant  Rankine 
source.  As  already  mentioned,  the  Rankine  source  method  customarily  adopts  the  upstream  finite-difference 
operator  to  satisfy  the  free  surface  condition  and  the  radiation  condition,  but  the  use  of  the  upstream 
finite-difference  operator  is  equivalent  to  introducing  an  artificial  viscosity.  This  viscosity  effect  attenuates  the 
wave  amplitudes  significantly  as  they  propagate  and  the  results  get  progressively  poorer  far  downstream.  There 
should  be  a  strong  effort  to  develop  a  computational  method  that  combines  the  excellent  near-field  capability  of 
the  Rankine  source  method  with  the  superior  far-fieid  behavior  of  the  Havelock  source  method. 

CONCLUSIONS 

A  higher-order  panel  method  for  the  linearized  free-surface  flow  has  been  developed.  The  capability  of  the 
present  higher-order  curved  panel  method  is  summarized  as  follows: 

1 .  The  cost  of  curved  higher  order  panels  in  the  computation  of  the  influence  coefficients  is  not  notably 

greater  than  the  cost  of  flat  constant  source  panels,  because  of  the  use  of  recursive  relations  from 

the  flat  constant  case. 

2.  The  results  from  the  use  of  curved  panel,  linear  source  terms  on  the  ship  hull  are  robust  with  respect 

to  the  panel  size  and  geometry. 

3.  The  convergence  rate  for  the  double  body  case  is  O(h^).  while  the  zerc'-order  flat  panel  method  gives 

0(h). 
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4.  The  free  surface  calculations  have  been  shown  to  be  robust  with  respect  to  changes  in  the 
longitudinal  and  lateral  free  surface  computational  domain  and  with  respect  to  the  panel  density. 

5.  New  partitioned  matrix  techniques  have  been  used  to  minimize  storage  requirements  and  allow  the 

number  of  body  and  free  surface  panels  to  be  up  to  2800  panels. 

6.  Extensive  use  of  advanced  vectorization  techniques  has  made  the  code  highly  efficient. 

7.  A  new  transom  stern  condition  has  been  developed  and  implemented  which  gives  excellent  results  for 

the  tested  models. 

8.  A  number  of  tests  with  the  Wigley  hull.  Series  60,  ATHENA,  and  a  large  combatant  hull  have  given 

excellent  comparisons  with  measurements  of  wave  profiles. 

9.  A  new  user-friendly  preprocessor  program  has  been  developed  by  Mr.  Steven  Fisher  of  DTRC  to 

provide  input  to  SWIFT. 

10.  Two  new  postprocessors  have  been  developed  which  give  graphical  results  of  both  th,;  wave  height 

on  the  ship  body  and  the  contours  of  the  wave  heights  over  the  computational  region  (see 
* 

Robinson,  for  example) 


RECOMMENDATIONS 

1.  Since  the  wave-making  resistance  is  a  small  quantity  whose  accuracy  depends  critically  on  the 
corresponding  accuracy  of  numerical  discretization,  an  approach  other  than  a  direct  numerical 
integration  of  the  pressure  should  be  considered. 

2.  In  order  to  fully  exploit  the  superior  potential  of  the  present  higher-order  method  over  the  zero-order 

approach,  the  free  surface  condition  should  be  recast  in  a  form  which  minimizes  the  artificial 
numerical  dissipation  and/or  numerical  dispersion. 


Robins,  C.R.,  "Ship  Wave  Color  Graphics  Using  01-3000."  Computer  Science  Senior  Project.  U.  of  North 
Carolina  at  Charlotte,  1986. 
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3.  An  accurate  prediction  of  near-field  flow  around  a  surface  ship  is  essential  for  the  prediction  of  wave 
breaking  near  the  bow  and/or  stern,  for  investigation  of  pressure  distributions  (possibly  cavitation 
inception)  on  a  sonar  dome,  for  tracing  bubble  sweep  along  the  hull,  and  for  the  alignment  of 
appendages.  Since  the  local  and  nonlinear  effects  close  to  the  hull  surface  are  important,  a  nonlinear 
formulation  of  the  free  surface  flow  is  needed  for  the  accurate  analysis.  Commonly  used  surface 
singularity  techniques  for  nonlinear  free  surface  problems  are  based  on  the  iteration  of  the  unknown 
surface  wave  topology,  in  which  the  adjustment  of  the  free  surface  shape  obtained  from  the  previous 
iteration  is  a  crucial  factor  in  the  convergence  and  stability  of  the  iterative  solution.  The  present 
panel  method  allows  the  use  of  higher  order  curved  panels  which  can  closely  approximate  the  actual 
surface  topology.  Further,  the  linear  variation  of  surface  singularity  in  the  higher-order  panel 
method  enables  self-induction  effects  to  be  included  for  each  panel.  The  higher-order  panel  method 
is,  therefore,  expected  to  make  the  iteration  solution  more  stable  than  that  of  existing  lower  order 
panel  approaches.  The  nonlinear  extension  of  the  present  numerical  scheme  will  be  the  subject  of 
separate,  ongoing  work. 
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APPENDIX  A 


COMPUTATION  OF  VELOCITY  POTENTIAL  i 

In  this  appendix  we  will  evaluate  the  integrals  of  Eqs.  (21).  (22).  and  (25).  With  the  relations  of  Eqs.  (15). 
(17).  (18),  and  (19),  the  induced  velocity  potential  at  P(x,y.z)  due  to  the  singularity  distribution  on  the  panel  L 
can  be  expanded  as 

♦  =  xl^(l,l)+l^(2.1)  1  +  ff^[  yl^(1.1)+l^(1.2)  ]  (Al) 

and 

♦  =  T^[xi.^(Li)+y2,i)]  +  y  yyi.i)+yi.2)  ] 

+  (1/2)t^^[  x^yi,l)+2xy2,I)+y3.1)  ]  +  (I/2)T^^[  y\(I.I)+2yyi.2)+yi.3)  ] 
+  xyyi,I)  +  yy2.1)  +  xyi.2)  +  y2,2)  ].  (A2) 

Here 


I^(M,N)  =  (-IMr)  {  H(M,N.1)  +  a[hH(M+2.N.3)  +  2xhH(M+ I  .N,3)] 

+  b[hH(M.N  +  2.3)  +  2yhH(M.N+ 1.3)] 

+  cIhH(M.N.3)]  }  (A3) 

and 

yM.N)  =  (l/4n)  {  hH(M.N.3)  +  aIH(M  +  2.N.3)  +  3h^H(M  +  2.N.5)+6xh^H(M+ 1  .N.5)] 
+  b[H(M.N + 2.3)  +  3h^H(M.N  +  2.5)  +6yh^H(M  .N  +  1 .5)] 

+  c[-H(M,N,3)+3h^H(M.N.5)]  }  (A4) 

where 


c 


2  ^  K  2 

ax  +  by  -  z 


(A5) 


(A6) 
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ax 

o 


+  by 
■^o 


H(M.N.K)  = 


f  f 


j  j 


(^x)^''(n-y) 

ic 


N-l 

-dwn 


p 


.i 


(^x)^  +  (n-y)^  +  h^ 


and  (’^q'Vq'O)  's  the  point  on  £  closest  to  (x,y,0)  (see  Fig,  2).  Since 

(3/3x)H(M.N.K)  =  -(M-1)H(M-1.N,K)  +  KH(M  +  1  .N.K+2) 
(3/3y)H(M,N,K)  =  -(N-l)H(M.N-l.K)  +  KH(M.N+ 1  .K  +  2) 
(3/3z)H(M,N,K)  =  -KhH(M,N.K+2)  , 


the  induced  velocities  are  obtained  as 

V*  +  d^[5j^(i.i)+y2.i)]  +  o^[yyi.i)+yi.2)  ] 

where 

T  (M,N)  =  [  J  (M.N)  ,  J  CM,N)  ,  J  (M.N)  1 
a  ox  oy  oz 

and 


(A7) 

(A8) 

(A9) 


(A  10) 


(All) 


(A  12) 


J  (M,N)  *  (-l/4n){  H(M+1,N.3)  +  3ah[H(M  +  3.N.5)  +  2xH(M  +  2.N.5)] 
ox 

+  3bh[H(M  +  l,N  +  2.5)  +  2y(M+l.N  +  L5)] 

+  3chfH(M+l.N.5))  I  (A13) 

J  (M.N)  -  (-l/4n){  H(M.N  +  1.3)  +  3ahtH(M  +  2.N+ 1 .5)  +  2xH(M+ 1  .N+ 1 .5)] 
oy 

+  3bh[H(M  .N  +  3.5)  +  2yhH(M  .N  +  2 .5)] 

+  3ch[H(M.N+1.5))  !  (A 1 4) 
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J^(M.N)  =  (-l/4n){  -hH(M.N.3)  +  afH(M  +  2.N.3)-3h  H(M  +  2.N.5) 

+  2xH(M+ I.N.3)-6xli  HfM+I.N.5)] 

+  b[H(M.N  +  2.3)-3h  H(M.N  +  2.5) 

+  2yH(M.N+1.3)-6yh  H(M.N+1.5)] 

+  c[H(M.N.3)-3h  H(M.N.5)]  }.  (A15) 

The  corresponding  expressions  for  the  induced  velocities  due  to  doublet  can  be  obtained  easily  by  diff¬ 
erentiating  Eq.  (A2)  and  then  using  the  relations  of  Eq.  (AlO). 

If  the  field  point  P(x,y,z)  is  located  at  a  large  distance  from  the  panel  surface  S,  the  expressions  (A3). 
(A4),  (A14),  and  (AIS)  can  be  further  simplified  using  a  series  expansion  with  respect  to  origin  of  the  local 
coordinate  system  For  large  values  of  |P|  compared  to  |Q|.  we  obtain 

I^(M.N)  =  (-1/411)  {  El/P  +  (E2  P)/P^  -I-  [-E3/2-t-3(P  E4P)/2]/P^  }  (A16) 

where 


El  =  C(M,N)  (AI7) 

E2  *  [  C(M-fI,N)  ,  C(M.N-H)  .  aC(M+2.N)-t-bC(M.N-l-2)  ]  (A18) 


E3  *  C(M+2.N)  C(M,N-b2) 


(AI9) 


E4 


fC(M-l-2.N)  C(M-t-l.N-l-l)  aC(M-t-3.N)  +  bC(M-t-l.N  +  2) 

C(M-I-1,N-H)  C(M.N  +  2)  aC(M-(-2.N-H)-l-bC(M,N-t-3) 

aC(M-b3.N)-H)C(M+l.N  +  2)  aC(M  +  2.N-)-l)+  bC(M.N-»-3)  0 

(A20) 


C(M,N)  >=  H(M,N,0)  evaluated  at  x  =  y  =  z  =  0.  (A2I) 

The  leading  terms  in  Eqs.  in  (A3).  (A4).  and  (All)  correspond  to  a  flat  panel  The  remaining  terms  having 
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coefficients  of  a,  b.  and  c  constitute  panel  curvature  effects.  The  corresponding  integrals  associated  with 
doublet  distribution  can  be  expressed  in  terms  of  H  integrals  with  different  combination  of  values  M.  N.  and 
K.  Furthermore,  we  can  show  that  there  exists  a  recurrence  relationships  between  H(M.N.K)  such  that  the 
computation  of  velocity  potential  and  its  derivatives  involves  only  simple  algebraic  manipulations.  The  detailed 
evaluation  of  the  H  integral  is  illustrated  by  Johnson^  and  the  results  are  summarized  in  Appendix  B. 

The  total  induced  velocity  potential  at  P(x..y..z.)  is.  therefore,  obtained  by  collecting  individual  con- 
-tributions  due  to  all  the  surface  panels  N. 


N 

E  ♦ 
j=l  •’ 


(A22) 


Using  the  the  relation  for  the  source  distribution,  Eq.  (20).  the  total  induced  velocity  potential  and  induced 
velocities  are  expressed  in  terms  of  N  singularity  parameters 


N 

j:  (€♦),, ff 
j=i 


u  J 


N 

V.  =  E  [  CV  .  CV  ,  CV  ]..  cr. 

'  i_i  *  y  z  U  J 


(A23) 


(A24) 


where  (C+)jj  and  [CV^,CV^,CV^].j  denote  the  influence  kernels  associated  with  the  j-th  singularity  parameter 
evaluated  at  the  i-th  control  point. 
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APPENDIX  B 
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Similarly  we  can  devise  the  reciircion  relations  for  the  computation  of  F  infeprals.  From  the  definition  of  F 
integrals  we  have  the  identities 

F(M  +  2,N.K)  +  F(M,N  +  2,K)  +  h^F(M.N,K)  =  F(M.N.K-2)  (B6) 

and 

Vj.F(M+l,N,K)  +  \)^F(M.N  +  1,K)  =  -aF(M,N.K).  (B7) 

Integration  by  parts  also  yields 

(K-2)\)^F(M,N  +  1.K)  -  (K-2)\)^F(M+1,N.K) 

=  Nv^F(M,N-l,K-2)  -  Mv^F(M-l.N.K-2)  +  E(M.N.K-2)  (B8) 

where 

E(M.N.K)  =  [  (e-x)^'’(tvy)^‘’/p*^  ].  (B9) 

The  recursion  relations  for  H  and  F  have  been  recombined  in  the  present  procedure  to  compute  needed  H 
and  F  integrals  efficiently  and  accurately.  For  example,  since  S  is  a  singular  surface  for  H.  the  reverse  of  the 
relation  Eq.  (B3)  has  been  used  for  the  field  point  in  the  interior  of  S.  Similar  care  has  been  exersised  for  the 
evaluation  of  F  integrals.  By  properly  accounting  for  the  singular  behavior  of  H.  it  was  possible  to  compute 
accurate  velocities  even  when  the  field  point  is  close  to  the  edge  or  corner  of  S. 
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I 


Fig.  1.  Coordinate  system. 


Fig.  3.  3*-D  perspective  view  of  the  surface  panel  network. 
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Rg.  5.  Wigley  hull  form  (from  Shearer  aruf  Cross,  1 965). 
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Relative  Errors  in  C„  Relative  Errors  in 


1.0  1.5  2.0  2.5  3.0  3.5 

Forward  and  Aft  Truncation  Boundary  (Xt/(L/2)) 


Fig.  6.  Convergence  test  for  X-direcHon  truncation  boundary. 
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Side  Truncation  Boundary  (Ys/(L/2)) 


Rg.  7.  Convergence  test  for  Y-direction  truncation  boundary. 
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Pig.  9.  niustration  of  body  and  free  surface  panel  network. 
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Absolute  Velocity 


Regular  Paneling  Random  Paneling 
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Wave  profile  (t/(UV2g))  Wave  profile  (t/(U  /2g)) 
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Pig.  11.  Wave  profile  comparison  for  Wigley  hull  at  Fn-0.250. 
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ng.  12.  Wave  profile  comparison  for  Wigley  hull  at  Fn-0.267. 
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Wave  profile  (C/(U  /2g)) 
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Pig.  13.  Wave  profile  comparison  for  Wigley  hull  at  Fn>0.289. 
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ng.  14.  Wave  profile  comparison  for  Wigley  hull  at  Fn-0.316. 
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Wave  profile  (t/(U  /(2g))  Have  profile  (i;/(U^/(2g)) 
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Rg.  15.  Wave  profile  comparison  for  Wigley  hull  at  Fn-0.354. 
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Rg.  16.  Wave  profile  comparison  for  Wigley  hull  at  Fn«0.408. 
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Rg.  17.  Wave-making  resistance  comparison  for  Wigley  hull. 


Fig.  18.  S-D  perspective  view  of  Wigiey  hull  generated  waves  at  Fn-0.408. 
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Hg.  22.  Wave  profile  comparison  for  Series  60  at  Fn«0.250. 
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Wave 


Longitudinal  axis  (2x/L) 


Rg.  23.  Wave  profile  comparison  for  Series  60  at  Fn>0.280. 
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Rg.  24.  Wave  profile  comparison  for  Series  60  at  F„»0.300. 
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Wave  profile  (C/(u  /(2g)))  Wave  profile  (c/(U  /(2g))) 
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Rg.  26.  Wave  profile  comparison  for  Series  60  at  Fn^O-SSO. 
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Rg.  27.  Wave-making  resistance  comparison  for  Series  60. 


Fig.  28.  Longitudinal  wave-cut  comparison  for  Series 
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LongltudlMl  axis  (2x/L) 

Rg.  32.  Wave  profile  comparison  for  Athena  hull  at  Fn-0.412. 


Wave  profile  (C/(D^/(2g)))  Wave  profile  (C/(U  /(2g))) 
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Hg.  34.  Wave  profile  comparison  for  Athena  hull  at  Fn*0.651 . 
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Fig.  35.  Wave  profile  comparison  behind  stem  at  Fn-0.484. 
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MODEL  A 
MODEL  X 


Rg.  36.  Body  plans  of  Models  X  and  A. 
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Fig.  37.  Body  plans  of  Models  A  and  C. 
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Wove  Height  Relotive  to  Boseline  Cft)  Wove  Height  Reiolive  to  Baseline  (ft) 
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Rg.  38.  Wave  profile  aiortg  the  hull  at  Pn-0.179. 
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Rg.  30.  Wave  profile  along  the  hull  at  Fn«0.238. 


54 


Wove  Height  Relotive  to  0oseline(ft)  Wove  Height  Relotive  to  Boseline  (ft) 
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Rg.  40.  Wave  profile  along  the  hull  at  Fn-0.298. 
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Rg.  41.  Wave  profile  alortg  the  hull  at  FnaO.179  for  Models  A  and  C. 
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Wove  Height  Relative  to  Baseline  (ft)  Wave  Height  Relative  to  Baseline  (ft) 
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Hg.  42.  Wave  profile  along  the  hull  at  Fn-0.238  for  Models  A  and  C. 
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Fig.  43.  Wave  profile  along  the  hull  at  Fn-0.298  for  Models  A  and  C. 
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FIq.  45.  Longitudinal  wave  cut  at  Fni4.238  (Probe  1). 
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Wove  Height  Relotive  to  Colm  Woler  (ft)  Wove  Height  Relolive  to  Colm  Woter  (ft) 
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Fig.  46.  Ljongitudinal  wave  cut  at  Fn«0J238  (Probe  2). 
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Fig.  47.  Lxxigitudinal  wave  cut  at  Fn«0238  (Probe  3). 
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Wove  Height  Relotive  to  Calm  Water  (ft)  Wave  Height  Relative  to  Colm  Water  (ft) 


Fig.  48.  Longitudinal  wave  cut  at  Fn-0.238  (Probe  4). 
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Fig.  40.  Longitudinal  wave  cut  at  Fn>0.238  (Probe  5). 


16  18  20  22 


59 


Table  1.  Comparisons  between  flat  panel  and  curved  panel  results. 


Relative  Error(^) 

No.  of 

Panels 

Flat  Panel. Constant  Source 

Curved  Panel. Linear  Source 

* 

X 

6 

y 

♦ 

z 

m 

B 

B 

2 

50 

3.43 

4.52 

1.85 

1.85 

0  64 

0.9tJ 

A  91 

0.91 

200 

I.7J 

1.25 

0.32 

0.32 

0  12 

0.20 

0.27 

0.27 

520 

M9 

0.66 

0.22 

0.22 

0.07 

0.17 

0.15 

0.15 

Note  :  Since  (he  body  is  axisymmetric.  one*<)uarter  of  body  panels  is  provided. 

The  relative  error  is  defined  as 

compMtrd  -  annlvtic 

relative  error  -  ABS( .  }  \  1 00 

anafxiic 


Table  2.  Tables  of  offsets,  Series  60,  Cb-0.60. 
(from  Todd,  1953) 

Halt  brMdthi  of  watarllna  givan  aa  traction  et  maximum  baam  on  oach  watarilna 


Modal  •  4210W 

W.L.  1.00  la  tha  daalgnad  load  watarllna 


rorabody  bilamatlc  coattlelant  •  O.HI 
Atiarbedy  pricmatlc  eoafticlant  •  0.M( 
Total  prlamatlc  coatllciont  •  O.IM 

Araa  aa 

fraction 

' Waiofilnoa  —  —  ■  . . . .  i.— i.,  et  max. 


>ta 

Tan. 

fP 

0.000 

M 

0.000 

1 

0.013 

1H 

0.011 

2 

0.024 

1 

0.0H 

4 

0.1M 

1 

OJTf 

« 

0.40 

7 

0.004 

1 

0.031 

f 

o.»« 

to 

1«I0 

11 

o.tts 

12 

0.M2 

11 

0.707 

14 

4.<Z2 

11 

0413 

It 

0.300 

17 

0.111 

It 

0.016 

1t» 

nos 

If 

O.CM 

IfVI 

0.010 

AP 

0.000 

Max  htH 

baam* 

0.710 

0.07t 

0.2S 

0.000 

4.000 

0.032 

0.OU 

0.0M 

0.0t2 

0.0K 

0.120 

0.127 

0.170 

O.lOf 

0.2M 

0.314 

0.410 

0.404 

0.M 

o.no 

0.723 

0.771 

OJM 

out 

0.036 

O.tM 

O.tTt 

1.000 

1.000 

0.M2 

0.000 

0.t22 

om 

0J2t 

0.012 

0.701 

0.701 

OfOO 

0.030 

0.413 

0.403 

92n 

O.IK 

0.103 

O.W 

0.120 

0.CH 

0.070 

0.020 

0.000 

0.000 

0.000 

0.iH 

0.106 

0.10 

0.000 

0.000 

0.041 

o.oc 

0.007 

0.000 

0.141 

0.100 

0104 

0.213 

0.l4t 

0J« 

OJOO 

0J30 

0.000 

0.011 

0002 

o.n4 

o.tot 

0.017 

0.071 

0.077 

O.0M 

1.000 

1.000 

1.BOO 

1.000 

1.000 

0.0M 

1.000 

0.|t2 

0.007 

0104 

0.0«3 

0.714 

0.007 

0.102 

0.730 

0.413 

0J41 

0130 

0J21 

0.1it 

0.210 

0.006 

O.IM 

0022 

0.033 

0.000 

0.000 

1.000 

1.000 

1.00 

116 

0.000 

0.020 

0.051 

0.070 

0.1CB 

0.133 

0.110 

0.166 

0126 

0170 

0.361 

0.440 

0.662 

0.007 

0.716 

0.764 

0.641 

0.162 

0126 

0.636 

6.676 

O.ltl 

1.600 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.004 

0.007 

0.176 

0.600 

0.637 

0.177 

0.667 

0.6X3 

4.736 

0.644 

0.636 

0.701 

0.436 

0.620 

0106 

0.630 

0.113 

0.410 

e.OK 

0170 

1.000 

1.000 

110 

araa  ta 

1.00  W.L. 

0.042 

0.000 

0.120 

O.Ott 

0.160 

0.066 

0170 

0.1V 

OHO 

0.112 

0.631 

0123 

0163 

0.676 

OIW 

0.630 

OM 

0.771 

0.040 

OM 

0.602 

0.M 

1.000 

0.660 

1.000 

1.000 

1.000 

0.M 

1.000 

o.tn 

IMO 

0.M 

0.666 

0.163 

0.164 

0.760 

0.676 

0.M 

0124 

0.446 

0.614 

OM 

0.760 

0.167 

O.OH 

0.166 

4.676 

0.040 

0.410 

0.004 

1.000 

*Aa  tracclen  at  maximum  lead  tyaiarlina  beam. 
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Table  3.  Offsets  for  tfie  high-speed  hull,  Athena. 


Sta. 

Tan. 

0.125 

03S 

0.50 

0  75 

1.00 

1.25 

1.50 

FP 

OSKtOO 

0.0000 

0.0000 

0.0000 

0.0000 

0.0048 

0.0185 

0.0347 

t/2 

0.0000 

0.0000 

0.0246 

03359 

0.0451 

0.0758 

a0989 

1 

0.0000 

0X1000 

0.0525 

0.0818 

0.0959 

0.1110 

0.1354 

0.1637 

1  1/2 

0.0000 

OXXKM 

0.0838 

0.1292 

0.1462 

0.1675 

0.1945 

03267 

2 

0.0000 

0.0000 

0.1162 

0.1766 

03035 

03257 

03S42 

0.2886 

3 

0.0000 

0.0377 

0.1955 

03813 

03104 

03398 

03711 

0.4081 

4 

OiMXK) 

0.1029 

03849 

03891 

04218 

0.4478 

0.4761 

05129 

5 

0i)000 

0.1972 

03989 

04992 

05280 

05643 

05776 

0.6078 

6 

OMOO 

0.3036 

0.4972 

0.6009 

0.6246 

0.6462 

0^00 

0.6990 

7 

OJOOOO 

0.4305 

03100 

03934 

0.7070 

0.7263 

0.7476 

0.7703 

8 

0.0000 

0J5918 

0.7262 

0.7783 

0.7830 

0.7967 

03156 

03345 

9 

OJOOOO 

0.7410 

0J8346 

03517 

0.8448 

0.8568 

03807 

0.8889 

fO 

IjOOOO 

OJOOOO 

0.9240 

03136 

03002 

0.9065 

03177 

0.9199 

n 

1.0000 

IjOOOO 

1,0000 

03671 

03420 

0.9381 

0.9457 

0.9488 

12 

IjOOOO 

0.83S3 

03619 

13000 

03762 

03660 

03684 

03699 

13 

OJOOOO 

0.4580 

03424 

13000 

03942 

03872 

03875 

0.9888 

14 

0.0000 

OjOOOO 

0376S 

03801 

1.0000 

1.0000 

1.0000 

1.0000 

IS 

OXMOO 

0.0000 

0.0581 

03113 

03865 

03939 

03946 

03953 

16 

OjOOOO 

0.0000 

03000 

0.7645 

03575 

0.9751 

0.8791 

0.9799 

17 

OjOOOO 

0.0000 

0.0000 

04870 

03227 

0.9478 

03505 

03517 

18 

OjOOOO 

0.0000 

03000 

03871 

03731 

03108 

03147 

03134 

18  1/2 

OjOOOO 

0.0000 

OjOOOO 

03000 

03545 

03926 

03920 

0.8899 

19 

OjOOOO 

OilOOO 

0.0000 

0.0000 

03345 

0.8695 

03711 

0.8669 

19  1/2 

0i>000 

0.0000 

03000 

03000 

03068 

03477 

03461 

0.8416 

20  <AP) 

OjOOOO 

0.0000 

0.0000 

03000 

0.8023 

03289 

03216 

03168 

Max.tulf 

beam 

0.0073 

03638 

0J6431 

0.7937 

03424 

1.0000 

1.0170 

13303 

Table  4.  Principal  dimensions  and  form  coefficients  of  Models  X,  A,  and  C. 


Model  X 


nii(Cti**t  aiHCNsioMS 

MoapiMCNSiONAi.  cncprtCKHTS 

UPfCTN 

900.00  ft  (27*. 11  «> 

0.9*9 

Cvp 

0.771 

0.009 

UHCTM  (WU 

9I0.«S  fc  (277.97  •) 

Cp 

0.979 

0.979 

L|/VW) 

0.79) 

•CAM  (•() 

1*1.10  ft  («).0I  •) 

Cwf 

0.9** 

0.*9* 

O.OOft 

eiAn  (Tt) 

99. M  ((  (10. fe  •} 

Cm 

0.«I9 

Cs 

2.**9 

U/Lwl 

0.212 

ntK 

0.00  (C  (♦  tow) 

0.97* 

( 

0.0*9 

Fl/Ul 

0.90* 

oisPucimiTT 

709M  1  (TIMS  t) 

0.919 

Uf9/«x 

*.*9* 

IT/Ur) 

0.992 

verrKO  wtfACi 

1279*2  ft  (iim.2  to  •> 

Cl 

0.999 

9,/T, 

9.990 

c» 

0.0(1927 

tc 

7.9*9 

0.709 

0.000 

s 

7.007 

u 

10.21*  eigr««« 

Ctfto 

0.992 

•t'9j 

0.000 

®  . 

*.799 

in 

20.970  togrtot 

CtfPA 

0.997 

Tt'7, 

0.000 

«/(.0IU)) 

99. *2 

Model  A 


putMClMi  niMCMinat 

aoeotNCitsio«AL  mpricicirrs 

LIHCTH  (Lfp) 

900.00  (c  (27*. 92  ■) 

Cg 

•  0.999 

C»p 

0.91* 

Ag/Al 

O.OOl 

UNCTH  (Utl) 

919.99  ft  (279.99  •) 

Cr 

•  0.992 

C»rp 

0.997 

l<7Uri 

•CAM  (li) 

12*. 21  t*  (99.47  •) 

Cpp 

•  0.9*9 

CVPA 

0.729 

Lg/Wl 

0.000 

DUPT  (T.) 

99. 91  fc  (10.92  «) 

CpA 

•  0.*l9 

Cg 

2.719 

U/kwl 

0.219 

niM 

0.00  ft  (*  tow) 

CpC 

-  0.979 

( 

0.0** 

ft/Lwl 

0.9(7) 

ei5PuccHP.trr 

70090  L  (7I2M  t) 

CfK 

•  0.99* 

Ul/li 

7. *21 

fT/Ul 

0.999 

VCTTKO  gUIPACe 

1292*7  M  f*  (119*9.9  to  •> 

Cl 

•  1.020 

•iTTj 

9.922 

c* 

0.00921 

7.297  togrt«« 

Cwr 

•  0.729 

%/*! 

0.000 

7.079 

9.929  togr««« 

•  0.997 

1.771 

20.009  togrto* 

CwTA 

•  0.929 

Model  C 
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